Abstract. We present a simplified model for linearized 
perturbations in a fluid with both differential rotation 
and differential vorticity. Without the latter the model 
reduces to the classical Shearing Sheet used in the descrip- 
tion of spiral density waves in astrophysical disks. Without 
the former it reduces to the /3-plane approximation, used 
in the description of Rossby waves. Retaining both, our 
model allows one in general to discuss the coupling be- 
tween density waves and Rossby waves, resulting in what 
- - - is known as the "corotation resonance" for density waves. 
\^ . Here we will derive, as a first application of this model, the 
' properties of Rossby waves in a differentially rotating disk, 
i; We find that their propagation is quenched by differential 
I . ■ rotation: after a limited number of oscillations, a Rossby 
wave collapses to a singular vortex, as fluid elements are 
■ sheared apart by differential rotation. In a kcplcrian disk, 
, this number of oscillations is always lower than one. We 
also describe how, in a similar manner, a vortex is sheared 
in a very short time. 
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1. Introduction 

Spiral density waves are well known in differentially ro- 
tating disks; they can be propagated by pressure alone, 
producing the Papaloizou-Pringle instability (Papaloizou 
and Pringle, 1984); or they can rely on long-range forces: 
self gravity, as in galactic disks or Saturn's rings (see e.g. 
Binney and Tremaine, 1987, and references therein), or 
magnetic stresses (Tagger et ai, 1990; Tagger and Pellat, 
1999). Their nature is that of compressional waves, i.e. 
sound waves modified by differential rotation and long- 
range forces. They can be excited by various interactions 
(in particular tidal forces), or linearly unstable. 

On the other hand Rossby waves are widely studied 
in the context of planetary atmospheres, with solid rota- 
tion (associated with the planet) but differential vorticity, 
due to the latitudinal gradient. Their nature is essentially 
that of a torsional wave, propagating vorticity. Their ba- 
sic properties are derived in a convenient model, the "/?- 
plane approximation" , where one neglects the gradients of 
all equilibrium quantities (including geometrical ones) ex- 
cept vorticity: the physics is thus reduced to the minimal 
description needed to describe Rossby waves. 

Symmetrically, density waves are described in the 
"Shearing Sheet" model, where the sole gradient consid- 
ered is differential rotation, and in particular the gradient 
of vorticity is neglected. These models are very convenient 
since they allow one to consider the essential properties of 
the waves, including their amplification. Global models, 
on the other hand, give more complete results but their 
complexity in general precludes a detailed understanding 
of the basic physics involved. However the Shearing Sheet 
as well as the /3-plane approximation neglect quantities 
which may convey important physical information. In par- 
ticular, in a disk, the vorticitjfl is: 
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^ hereafter, without any risk of ambiguity, we will consider 
only what is actually the vertical component of vorticity 



where is the rotation frequency, and k the epicyclic fre- 
quency given by: 

= Afl^ + 2nn' r (2) 
where the prime notes the radial derivative. Thus in a kep- 
lerian disk (fi ~ r"'^/^) one has W = so that its gra- 
dient is comparable to differential rotation and should not 
a 'priori be neglected. Symmetrically, in planetary atmo- 
spheres, zonal circulation is often fast enough that differ- 
ential rotation might play an important role in the prop- 
agation of Rossby waves. 

In this paper we show how, in the cylindrical geometry 
of a disk, the derivation of the Shearing Sheet model can 
be made more exact, while retaining its essential simplic- 
ity. This results in a set of equations with two parame- 
ters, the gradients of rotation frequency and of vorticity. 
The model reduces to the /3-plane approximation or to 
the classical Shearing Sheet (CSS) respectively, when ei- 
ther parameter vanishes. In the former case one recovers 
Rossby waves, and in the latter case usual spiral density 
waves. In the general case where both gradients are re- 
tained, the model allows one to describe how these waves, 
loosing their identity because of differential rotation, be- 
come coupled: density waves generating Rossby waves as 
they propagate, and reciprocally Rossby waves spawning 
density waves as they are sheared by differential rotation. 
The corresponding exchange of energy and angular mo- 
mentum between the waves affects the linear stability of 
the disk, as discussed below. For this reason we call our 
model, appropriate to describe this exchange of energy 
and momentum, the Rossby Shearing Sheet. 

The details and consequences of this coupling will be 
presented in separate publications. Here we will discuss 
only the basic properties of the model, in section ^, and 
present in section ^ a first application: we will give an 
exact solution in the incompressible limit, describing the 
propagation of Rossby waves in presence of differential 
rotation. We find that after a finite number of oscilla- 
tions a Rossby wave always collapses to a singular vor- 
tex, as fluid elements participating in the wave motion 
are sheared apart by differential motion. We also discuss 
in section ^ how differential rotation shears a vortex in a 
very short time, as found in recent numerical simulations. 
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The last section will be devoted to a short discussion 
of these results, their consequences and their connection 
with already existing ones. 



2. The Rossby Shearing Sheet 

Our goal is to derive in a rigorous manner a form of the 
Shearing Sheet model, that would allow us to take into 
account the vorticity gradient, and thus the physics of 
Rossby waves. In practice this will result in a set of equa- 
tions going smoothly from the classical Shearing Sheet 
(CSS, no vorticity gradient) to the /3-plane approximation 
(no differential rotation). 

We thus start from the linearized Euler and continuity 
equations, written in cylindrical geometry: 
d 

— Vr - 2nVi} = fr (3) 

dt 
d_ 
dt 



disappear if we retained self-gravity) , while 



+ WVr = 

p = -pqV ■ V 



(4) 
(5) 



where f = F/po, po is the equilibrium density, and the 
force F is left unspecified; we have assumed, in the spirit 
of the CSS, a disk with constant density and temperature. 
We have defined: 
d d . ^ 

where m is the azimuthal wavenumber. 

We use logarithmic coordinates with 
s = Inr 
and define 



u 

V 

D 
R 



rVr 



— 7i + imv = r V • V 

OS 

d 

—V — imu = (V x V) 

OS 



Thus D and R represent the compressional and vortical 
components of the velocity field. 

Combining equations (0-B|) we get: 



d 

— D - 2nR - 2n'(v - imu) = r^V ■ f 
dt y J J 



dt^ 



WD + W'u = r^V X / 



(6) 
(7) 



d 



si"-—. (8) 

where the ' now defines the derivative with respect to s, 
and h = p/po- These equations are exact in cylindrical 
geometry. In MHD disks threaded by a vertical magnetic 
field (see Tagger & Pellat 1999) the term in V x / gives 
access to the physics of Alfven waves. Here we restrict 
ourselves to conventional fluid disks with only pressure 
stresses, so that this term disappears (it would also 



V/ 



f2 



where cs is the sound speed, and we have assumed for 
simplicity an isothermal equation of state. By writing the 
equations for D and R, we have explicitly displayed the 
term in W' in equation (|^), which does not appear in 
the usual derivation of the shearing sheet. Thus we can 
now use the same techniques as in the CSS, considering 
as constant all the coefficients in the left-hand sides except 
for the rotation frequency, linearized around the fiducial 
radius s = (which can conveniently be taken as the 
corotation radius rg, if one treats a monochromatic wave). 
Thus we write: 
d d 

di^dt^ *"i(f^o + 2Ams) (9) 
where Hq = fl{s = 0), and A = f2'/2 is Oort's constant. 
We perform a Fourier transform in s (note that then the 
radial wavenumber k is non-dimensional) , so that the term 
proportional to s gives a derivative in k. We work in the 
frame rotating at the frequency Qq, so that finally the time 
derivative becomes 
d d d 

We use the exact relations: 

—ikD + imR —imD — ikR 

where 

2 ; 2 , 2 

q = k + m 

and get a third order diff'erential system: 



dt 



d_ 

dt ' 
d_ 
dt 



2Am 



dk 



D 



2nR 



m D + (m^ - ik)R] + q'^c^h 



2Am 



2Am 



d_ 

dk 
d_ 
dk 



R 



h = 



-WD- 



"'I 



{kD ~ mR) 



(11) 
(12) 

(13) 
(14) 



where we have defined 

(note that a = ^ in a keplerian disk). In the CSS, the last 
term disappears from the right-hand side of equation (|l2|) , 
allowing us to reduce the system to second order with 



h = WriR 



(15) 



This describes the Papaloizou and Pringle (1984, 1985) 
instability; with additional forces it would describe self- 
gravity driven or magnetically driven spiral instabilities. 
Note that equation (|l^) is the linearized form of the 
more general conservation of specific vorticity (also called 
vortensity in this context): 
(V X V)/p 

but that a model involving equilibrium radial density gra- 
dients would be much less tractable than the present one; 
thus for simplicity we will stick (in the spirit of the CSS 
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and /3-plane approximation) to the description of a flow 
with no equiUbrium gradients besides rotation and vortic- 
ity. For comparison with more complete models however, 
we will need to keep in mind this difference. 

The third order system must usually be solved numer- 
ically: the coupling between the second-order system and 
the third equation, proportional to a, describes the ex- 
change of energy and angular momentum between the 
density and Rossby waves. However we can get inter- 
esting results already in the incompressible limit: setting 
h = D = 0, the system reduces to: 
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R = —im—R 



(16) 



Let us first consider solid rotation with differential vor- 
ticity (A = 0, a 7^ 0): equation ( p^ ) gives the dispersion 
relation of Rossby waves, in the /3-plane approximation: 

oj = ^ (17) 

Going back to the full compressible equations, and re- 
taining both A and a, we use the little-known, but power- 
ful formalism of Lin and Thurstans (1984): we note that 
t appears in equations (11- 1^) only through the opera- 
tor in the left-hand sides. We can thus separate variables, 
writing: 

$(fc, t) = ^{k)F{k - 2Amt) (18) 
where $ is the vector (_D, R, h). Then the function 
F{k — 2Amt) can be factored out and disappears from the 
equations; it is an envelope function, given by the initial 
conditions. Two choices of F have an important physical 
meaning: F = Cte corresponds to normal modes (stand- 
ing wave patterns). This can easily be seen by writing the 
inverse transform: 

$(s, t)^ / dk e'''-' $(fc) F{k - 2Amt) (19) 

so that, with F = Cte, $ does not depend on t, in the 
frame rotating at the frequency f2o- This corresponds to 
a standing wave pattern, i.e. a normal mode, for the dis- 
crete set of wave frequencies (corresponding to a choice 
of ro) such that the appropriate boundary conditions are 
fulfilled. 

In the same manner, one finds that F ~ d{k — 2amt — 
ko) corresponds to a single wave, launched with k — ko 
at t = 0. This allows one to study the transient fate of 
a single wave propagating in the disk, as differential ro- 
tation causes its radial wavenumber k (selected by the 
delta function in F) to evolve with time, so that the wave 
pattern changes from tightly wound leading (with k large 
and negative) to open (fc ~ 0) to tightly woimd trailing 
(fc +oo). 

For simplicity we will neglect the " below. Equations 
(0-u3) become, for any choice of F: 



d 

2Am—D 
ok 

\A 



2nR 



-2Am—R = -WD + i^(kD-mR) 
ok q-^ 

d 1 
~2Ani—h = — 



— [-(fc -I- i)mD + {m^ - ik)R] + q Cgh 



(20) 



(21) 

^dk-^ rr . "^'^ 

In equation (|2l|), the term proportional to a now describes 
an exchange of specific vorticity with the equilibrium fiow. 
With a = 0, the system can be reduced to a second or- 
der ODE, in which self-gravity and magnetic stresses are 
easily incorporated. This system has thus been the ba- 
sis to analyze the amplification of spiral waves driven by 
self-gravity, by magnetic stresses, or by pressure forces. 

As discussed above, on the other hand, in the in- 
compressible limit equation ( |2l] ) represents the propaga- 
tion of Rossby waves. In general the term proportional 
to a in equation ( ^ thus represents a source of Rossby 
waves, generated by the propagation of density waves. The 
Rossby wave exchanges energy and angular momentum 
with the spiral, thus causing amplification or damping. 

The dependence of this mechanism on the vorticity 
gradient, and the singularity of Rossby waves which will 
be discussed in the next section, allow us to identify it with 
the "corotation resonance", known to occur for density 
waves: this has been discussed in particular by Panatoni 
(1983) for self-gravity driven spirals, by Papaloizou and 
Pringle (1985), Papaloizou and Lin (1989) and Narayan 
et al. (1987) for the Papaloizou-Pringle instability, and 
more recently by Tagger and Pellat (1999) in the case of 
magnetically driven spirals. In the first two cases it was 
found to make only a minor difference on the properties 
and growth rate of the instabilities. 

In the latter case on the other hand it was found 
both quite efficient (as compared with the weak insta- 
bility found in the CSS) and promising: the spiral wave 
grows by extracting energy and angular momentum from 
the disk, and transferring them to the Rossby vortex at 
the corotation radius; the vortex twists the footpoints of 
magnetic field lines threading the disk, and this torsion 
will in turn be propagated as Alfven waves to the corona 
of the disk, where it might energize a jet or an outflow. 

As the CSS or the /3-plane approximation, our model 
can retain its essential simplicity only at the cost of dis- 
carding the gradients of equilibrium quantities such as the 
density or magnetic field. It is thus important to remember 
the more complete result obtained in the above-mentioned 
works : they have shown that the corotation resonance de- 
pends in fact on the gradient of W/ p (the specific vortic- 
ity) for waves driven by pressure or self-gravity, or W p/ B"^ 
(where is the equilibrium magnetic field) for waves 
driven by magnetic stresses. 

3. Rossby waves with differential rotation 

These applications of the coupling between spiral density 
waves and Rossby waves must be discussed in the spe- 
cific context associated with a given disk model, and we 
defer this to separate works. In particular we will make 
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Flg.JLf Real (solid) and imaginary (dots) parts of the solution, equa- 
tion (|23|) for very weak differential rotation (a/2 Am — 20). The finite 
imaginary part at — * /oo is responsible for the singular behavior of the 
wave as i — * oo. 



use of the fact that these equations can straightforwardly 
be extended to include vertical motions and the vertical 
gradient of equilibrium quantities, in order to study the 
vertical structure of the waves. 

Here, in order to illustrate the use of the Rossby Shear- 
ing Sheet, we limit ourselves to present an exact solution 
for the propagation of Rossby waves in a differentially ro- 
tating disk, in the incompressible limit. This result is new, 
although it might have been obtained much more simply 
from the conservation of specific vorticity. We believe that, 
given the simplicity of the system of equations ( pO| - ^ ), 
it can be applied straightforwardly to a planetary atmo- 
sphere with zonal circulation. 

In the incompressible limit, {D — h — 0) equation ( pT| ) 
has the exact solution: 

/Co exp 



R 



arctan 
2Am m 



(23) 



Let us first consider this solution with the envelope 
function: 

F{k - 2Amt) = Cte. (24) 
The inverse Fourier transform gives: 



R{s,t) = i?o 



dk exp 



ks 



a k 
— - — arctan — 
2Am m 



(25) 



which is independent of i (i.e. a normal mode solution if, 
by the choice of the frequency uj — mflo, we manage to 
verify the appropriate boundary conditions). 

We will first discuss a case with weak differential rota- 
tion, j4 <C a. Figure (|l]) shows the integrand, in equation 
(p^), for a/2 Am = 20. The exponent in equation ( psf ) 
varies rapidly when k m. The method of stationary 
phase tells us that at a given radius s the main contribu- 
tion to the integral in equation ( p5| ) comes from the value 
of k such that the exponent is stationary: 

-+ 2A(/+m^) ^» (2^) 
Using equation (^) we see that this gives a local 
wavenumber fc(s) such that it verifies the dispersion rela- 



tion of Rossby waves, Doppler-shifted by the equilibrium 
rotation: 



U! 



m 



ma 



(27) 



with UJ = mflo- The wave propagates between s — (coro- 

tation, at large k) and a turning point (fc = 0) at: 
a 

~2i;^ 

Let us now consider the solution with 
F = d{k- 2Amt) : 

the S function has the effect that the full solution R{k,t) 
explores the separated one, equation(p3|), while k varies 
from tightly wound leading (fc — > — oo) to tightly wound 
trailing (fc ^ oo) because of the shear. The inverse Fourier 
transform gives: 



R{s,t) = Ro 



dke''"'6{k - 2Amt) 



exp 



a 



k 

arctan — 
2Am m 



— Rq exp 



a k 
— - — arctan — 
2Am m 



(28) 



where now k — 2Amt. Again in the limit where differential 
rotation is weak, A ^ a, the stationary phase tells us that 
at a given s the perturbation is strong only when k has 
the value given in equation (p6|), and at the correspond- 
ing time: we now have a description of a travelling Rossby 
wave, propagating from tightly leading (fc large and neg- 
ative) at corotation at t = — oo, to the turning point (at 
t = 0) and back to corotation as tightly trailing (fc large 
and positive). 

A qualitatively new and important result is obtained 
by considering the evolution of the wave at corotation: as 
k varies the arctangent in equation ( p8| ) varies from —tt/2 
to 7r/2, so that the exponent varies between ±TTa/AAm. 

This means that at corotation the Rossby wave oscil- 
lates only 

rv 

(29) 



■illatii 



AAm 



times before it gets sheared away by differential rotation. 
H always evolves, as t oo, to a singular vortex (a jump 
in at s — 0). This final evolution will occur even for 
arbitrarily weak differential rotation. 

At s ^ the solution still oscillates with time because 
of the ks term in the exponent: this is only due to the 
fact that a feature which is stationary at s = is seen 
Doppler-shifted by differential rotation by an observer at 
a different radius. These results are still valid (since the 
solution, equation ( p3| ) is exact) when differential rotation 
is not weak, as in a keplerian disk, but their WKB inter- 
pretation in terms of propagating waves becomes blurred. 

At this stage we have thus described, in an equivalent 
(by the method of stationary phase) to a WKB model, 
the propagation of a single Rossby wave in a disk with 
differential rotation. 

In the normal mode solution {F = Cte), at s = 0, 
the initial and final values of the integrand in equation 
( psi) are different, also giving by inverse Fourier transform 
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Fig. 2. Evolution of a vortex in a differentially rotating disk. Contour 
plots, as well as cross sections at ^ — and s — 0, are shown as time 
varies from to tt, i.e. during a half rotation period. 



a singularity at corotation. In the standing pattern ob- 
tained with this choice of the Rossby wave becomes a 
stationary vortex. 

In general the envelope function F is given by the ini- 
tial conditions, so that if they are regular initially {e.g. 
describing a localized wave packet) the Fourier transform 
will be well-behaved at infinity. The evolution will never- 
theless lead to the same conclusion. 

In the more general case where A is not small, this 
evolution will be very rapid: in a keplerian disk {A = a) 
the wave will describe only a fraction l/4m of a cycle 
before it degenerates! This is the reason why Rossby waves 
cannot be seen in disks with a realistic vorticity gradient. 
It is also the reason why the corotation resonance is more 
efficient at low m, so that the CSS can be considered as a 
convenient approximation for m 3> 1. 

4. Evolution of vortices 

Recent numerical studies (Bracco et ai, 1998; Godon 
and Livio, 1999, 2000; Davis et ai, 2000) have addressed 
the fate of vortices in protoplanetary disks. This was 
prompted by the work of Barge and Sommeria (1995), 




Fig. 3. 

details. 



Same as figure 3, with color-coded contours showing more 



who had shown that a vortex would favor the forma- 
tion of planets, by capture and coagulation of dust par- 
ticles. These studies show that, as expected, vortices are 
rapidly sheared apart by differential rotation unless they 
are strong enough that non-linear effects allow them to 
survive for a significant time. This is obtained by Godon 
and Livio (1999, 2000) as embedded cyclonic and anticy- 
clonic vortices. Davis et al. obtain longer lived anticyclonic 
vortices. The difference probably lies in the initial condi- 
tions given. In both cases this requires such strong ampli- 
tudes (locally reversing the rotation velocity gradient in 



Tagger: Rossby waves and vortices 



7 



the work of Davis et al.) that the mechanism which could 
create such vortices (supersonic if they are larger than the 
disk thickness) remains mysterious. 

Here we will give a simple description of the linearized 
evolution of a vortex. We start from an initial gaussian 
vorticity perturbation, whose Fourier transform is: 
R°{m,k) = e-^'^'/"'o+kVkl) (30) 

with mo = ka — 7 for the numerical solution displayed. 
Using the previous results, we find that the envelope func- 
tion for each Fourier component is given by: 



F„,{k) = i?°(m,fc)exp 



a k 
i — - — arctan — 
2Am m 



(31) 



giving the time evolution of the vorticity (after a conve- 
nient change of variables, fc — > fc + 2Amt, in the inverse 
transform) : 



R 



gim(i?+2Ats)g- 



exp 



2Am 



, k 

arctan( — 
m 



2 At) ~ arctan(— ) 



(32) 



which is easily computed (as it would be with any other 
choice of initial condition). The result is displayed in fig- 
ures I and |. It shows that, as expected, the vortex is 
rapidly sheared away by differential rotation, reducing to 
the vorticity sheets described by Davis tt al. (2000). It 
applies (since this is linear theory) to positive or negative 
vorticity, i.e. cyclonic or anticyclonic vortices, and cor- 
responds to the evolution observed for all but the high- 
est amplitude vortices studied in non-linear simulations. 
An obvious extension would be to include compressibility, 
i.e. return to the full third-order system of equations, and 
solve it numerically for the equivalent of the incompress- 
ible solution, equation (^3|). 

5. Discussion 

We have derived, from the equations describing a differ- 
entially rotating disk, a model (which we call the Rossby 
Shearing Sheet) which retains the intrinsic simplicity of 
the conventional Shearing Sheet (obtained with only dif- 
ferential rotation) or the /3-plane approximation (obtained 
with only differential vorticity) while including both gra- 
dients. The model results in a third order system of ordi- 
nary differential equations, rather than second order as in 
the CSS or first order in the /3-plane approximation. This 
means that, besides inward- and outward- propagating spi- 
ral density waves, it incorporates the physics of Rossby 
waves. 

In a first application, we have used incompressibility to 
separate the propagation of Rossby waves; we have shown 
that, for any ratio Aj a between differential rotation and 
differential vorticity, the wave always collapses after a fi- 
nite time to a singular vortex (a jump in at its coro- 
tation radius). In a realistic disk, this evolution is so fast 
that it occurs before the wave has completed even one cy- 
cle of oscillation. In a planetary atmosphere, on the other 



hand, zonal circulation should result in a lower ratio Aja 
and a slower evolution of the wave, but still in its collapse 
to a singular vortex in a finite time. 

We have also described the linear evolution of a vortex 
in a keplerian disk, showing (as observed in non-linear 
simulations) that it reduces to a thin vorticity sheet in a 
very short time. This confirms that, unless a very strong 
non-linear mechanism can be found to maintain a vortex 
against shear, there is little hope of forming planets in this 
manner. 

Lifting the constraint of incompressibility, we find that 
the density waves and Rossby waves are no more inde- 
pendent: this means that they are coupled by differential 
rotation in the disk. Waves loose their identity (the so- 
lutions are no more separated), so that a density wave 
generates a Rossby wave, and reciprocally. This process, 
in which energy and angular momentum are exchanged 
leading to growth or damping of the waves, is known as 
corotation resonance (but was not analyzed in terms of 
coupling to the Rossby wave) in studies of the spiral insta- 
bility, with or without the long-range action of self-gravity 
or magnetic stresses. In the latter case it has recently been 
found to make an important contribution, leading to what 
we have called an Accretion-Ejection instability. In recent 
works, Lovelace et al. (1999) and Li et al. (2000) have 
also found an instability of unmagnetized disks, when the 
specific vorticity profile has an extremum. 

Actually our result already allows us to give a physi- 
cal interpretation for the sign of the corotation resonance: 
we found that a Rossby wave can propagate only between 
corotation and a radius Smax = —a/2Am^, i.e. inside coro- 
tation if a and A have the same sign: since a wave propa- 
gating inside corotation has negative energy and angular 
momentum (see e.g. Collett and Lynden-Bell, 1987) we 
conclude that a spiral wave propagating inside corotation 
will be damped by exciting a Rossby wave. On the other 
hand we have mentioned that in fact the relevant gradi- 
ent for Rossby waves is that of specific vorticity, W/T,, or 
of WT,/B^ in a magnetized disk. Thus if these gradients 
change sign, they will also change the sign of Smax so that 
the Rossby wave will now propagate beyond corotation, 
and will have a positive energy. This means that exciting 
it will amplify a spiral wave. 

Future work will be dedicated to using this model for 
a deeper analysis of the coupling between density and 
Rossby waves, in accretion disks and in planetary atmo- 
spheres. 
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